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1.0 Scope of Work 


This report describes the work performed by the University of North Dakota 
(UND) under NASA Grant NAG5-1 1509, titled “Airborne In Situ and Ground-based 
Polarimetric Radar Measurements of Tropical Convection in Support of CRYSTAL- 
FACE.” This work focused on the collection of data by two key platforms: the UND 
Citation II research aircraft and the NASA NPOL radar system. Prof. Michael Poellot led 
the deployment of the Citation and Dr. Paul Kucera led the operations of the NPOL radar 
in collaboration with John Gerlach (NPOL facility manager, NASA Wallops Flight 
Facility) and Rob Cifelli (Colorado State University). 

As stated in the NRA, the CRYSTAL-FACE (C-F) mission addresses several key 
issues from the NASA Earth System Enterprise, including the variability of water in the 
atmosphere, the forcing provided by tropical cirrus and the response of the Earth system 
to this forcing. In situ measurements and radar observations of tropical convection, cirrus 
clouds and their environment are core elements of C-F. One of the primary issues that C- 
F is addressing is the relationship of tropical cirrus anvils to precipitating deep 
convection. Because convective clouds act as a source region for tropical cirrus anvils, it 
is important to obtain observations of the microphysical and spatial properties of the 
parent convective cells through the entire life cycle of the cloud system. 

The in situ measurements from C-F are being used to validate remote sensing of 
Earth- Atmosphere properties, increase our knowledge of upper tropospheric water vapor 
and its distribution, and increase our knowledge of tropical cirrus cloud morphology and 
composition. Radar measurements, especially polarimetric diversity observations 
available from the NASA NPOL radar, are providing essential information about the 
initiation, modulation, and dissipation of convective cores and the generation of 
associated anvils in tropical convection. Specifically, NPOL radar measurements contain 
information about convective intensity and its vertical structure for comparison with 
thermodynamic and kinematic environmental measurements observed from soundings. 
Because of the polarimetric diversity of NPOL, statistics on bulk microphysical 
properties can be retrieved and compared to the other characteristics of convection and 
associated cirrus anvils. 

In summary, the central objectives of this proposal were to deploy the UND 
Citation research aircraft as an in situ sensing platform for this mission and to provide 
collaborative analyses of the recorded data and to deploy the NPOL radar to observe the 
characteristics of cirrus and parent convection. 

2.0 Objectives 

Citation 

The following are the proposed work objectives related to Citation operations: 

1. Integrate guest instruments on the Citation as able and as selected by the CRYSTAL- 

FACE science team to help meet program objectives; 


2. Operate the UND Citation research aircraft during the CRYSTAL-FACE 

observational period, including measurements of cloud microphysics, state 
parameters and winds; 

3. Process data recorded on the Citation DAS and provide the data to the program 

archive; 

4. Compile mission summaries for all flight missions; and 

5. Perform collaborative analyses of the Citation data with other CRYSTAL-FACE 

investigators. 

NPOL Radar 

The following are the proposed work objectives related to NPOL operations and 
data analysis: 

1 . Field phase of CRYSTAL-FACE - data collection; 

2. Quality control of NPOL radar data; 

3. Development of precipitation fields and polarimetrically tuned Z-R relationships from 

NPOL radar data; and 

4. Hydrometeor identification from polarimetric signatures examine the evolution of the 

convective systems observed by NPOL for priority cases 

3.0 Work Accomplished 

3.1. Citation 

3.1.1 Integration of guest and new instrumentation - CFDC, CVI, CIN, CPI, HVPS, TDL 
The suite of instrumentation to be flown on the Citation was determined through 
discussions with the CRYSTAL proposal review team. The final package (Table 1) was 
based on the measurement objectives and the capacity of the aircraft to carry and power 
the desired equipment. Some of these instruments were provided by guest scientists. The 
mission of the Citation to provide in situ cloud microphysical measurements was 
reflected in the choice of instruments. 

Preparation for this field deployment actually began in January, after the grant 
award. Integration of this suite was a complex task in terms of number and nature of the 
instruments, and three of them had never been flown on the Citation before (CFDC, CIN, 
TDL). This suite also represented an unprecedented impact on the weight, space and 
electrical capabilities of the Citation, requiring several modifications. These included 
design and construction of a new instrument rack, the addition of electrical wiring and 
inverters for research power, and the design and construction of two new inlets for 
sampling. All aircraft modifications also required engineering and FAA approval. The 
design of the rack and one of the inlets, along with the engineering and FAA approvals, 
were all accomplished by UND. A contractor performed rack construction and the other 
inlet was built by NCAR, where a modification to the mount of the CVI was also done. 
The wiring and inverter installation was done during scheduled aircraft maintenance in 
April. 


3.1.2 Deployment in Key West 

The Citation was ferried to Key West Naval Air Station on June 30 and returned 
to Grand Forks on July 3 1 . Research flight operations began on July 3 and a total of 13 
research missions were flown during the experimental period. A list of these flights and 
primary objectives is given in Table 2. Overall, this was a highly successful measurement 
mission for this platform. Instrument reliability was high; excellent coordination and 
support was provided through pre-flight briefings and real-time coordination by the crew 
at NPOL. As seen in Table 2, a variety of conditions were sampled according to the 
science plan of the experiment. 

3.1.3 Data Processing and Archival 

All flight data were processed after each mission upon landing. These processed 
data were then uploaded to the experiment archive. These were considered preliminary at 
that point, with quality control and final instrument calibration activities still in progress. 
After the field program the processing proceeded along several fronts in parallel with 
analysis efforts, including quality control and archival. 

The first step in quality control was to correct or flag known bad data associated 
with instrument problems noted during the flight missions. Instrument calibration issues 
were also addressed. Step two consisted of an extensive manual edit of the data. All 
recorded parameters were plotted and reviewed for reasonableness. Case study analyses 
also provided close scrutiny of select portions of the data. Measurements appearing 
abnormal or out of range have been flagged with missing data values. 

Archive formats for C-F data were established in a data exchange plan for this 
project (http://esDoarchive.nasa.gov/info/crvstalftdataex plan.htmD . The Citation data 
were grouped into four categories to meet plan requirements: meteorological data (MET), 
cloud physics data (CLD), navigation data (NAV) and other observations (MIS). An 
additional file containing derived particle size distributions from the PMS FSSP probe 
was also created for each flight (FSS). A full description of data processing techniques, 
file formats and related information is presented at 

http://aerosol.atmos.und.edu/citation/index.html . During the archival process, quick-look 
plots were generated for select parameters and may be viewed at 
http://aerosol.atmos.und.edu/citation/gallerv/Crvstal-FACE . All data have been 
submitted to the ESPO archive. 

3.1.4 Mission Summaries 

Additional quick-look information has been compiled into Flight Summaries for 
each of the Citation missions during C-F. Each flight was broken down into segments and 
summary information compiled for each segment. These summaries have been added to 
the quick-look plot gallery described above. The summaries are attached to this report as 
Appendix A. 

3.1.5 Collaborative Studies 

Several collaborative studies have been completed using Citation data from C-F 
to address key science questions. Of particular interest were observations taken at cold 


temperatures in a convective core on July 18 and those made in clear air on July 28 and 
July 29. The core penetration on July 1 8 was made at a temperature of approximately -35 
°C at an altitude of 10 km. Updrafts exceeding 20 m s' 1 were encountered along with 
regions of ice and supercooled liquid water droplets. These observations have led to the 
conclusion that in vigorous convection, heterogeneous ice nucleation processes are not 
sufficient to freeze all of the supercooled cloud droplets before homogenous freezing 
occurs (Heymsfield et al., 2004). In a separate analysis of the July 18 observations, 
environmental aerosol concentrations and model simulations point toward a significant 
role of mid-tropospheric aerosols in ice formation. Fridland et al. (2004) attribute the 
high concentration of small ice crystals in this cloud to the entrainment of aerosols which 
activate and later freeze homogeneously. This has implications for the possible effects of 
long-range transport of pollutants. Both studies considered the importance of 
homogeneous freezing to the radiative properties of anvil clouds. 

The flights of July 28 and 29 encountered a layer of Saharan dust. Measurements 
with the CFDC by DeMott found the dust particles to act as ice nuclei, with high 
concentrations of >1 cm' 3 . (DeMott et al., 2003). Observations with ground-based lidar 
confirmed the presence of the dust layer. Depolarization measurements of an altocumulus 
cloud layer indicated that the dust acted to generate ice crystals within (seed) the 
altocumulus at relatively warm temperatures (Sassen et al., 2004). Based on these 
observations it appears that the Saharan dust has similar properties and ice nucleating 
characteristics as Asian dust. This suggests possible large-scale climatic influences of 
dust with implications for changes due to anthropogenically-induced desertification. 

Joint Aircraft/Radar Study 

A study of the microphysical and radar properties of clouds occurring on July 16 
is underway. Preliminary results were presented at the International Conference on 
Clouds and Precipitation in Bologna, Italy (Poellot et al. 2004). 

3.2 NPOL Radar 

3.2.1 Field Operations in South Florida 

The NPOL radar was deployed in Southwest Florida about 10 km East of the 
Everglade City field site (Ochopee, Florida) near the of June. NPOL was operational on 
July 1 with the first day of 24 hr operations occurring on July 5. NPOL completed 
operations on July 29. Operational support was split between John Gerlach and Pi’s 
Kucera and Cifelli, with each PI in the field for approximately 3 weeks. Graduate and 
undergraduate students were invaluable in assisting in radar operations and preliminary 
analysis in the field. The radar operated almost continuously during this time supporting 
aircraft operations and collecting valuable information on the temporal evolution and 
three-dimensional structure of the convection observed during C-F. NPOL operated in a 
full- volume scan mode during the daytime hours and/or during aircraft missions and in a 
3-tilt rain scan mode during all other periods. NPOL operated every day during the C-F 
field campaign, but had some minor antenna problems that mainly occurred during July 
10-12. Intermittent air conditioning problems during the latter half of the campaign 
resulted in a reduced scanning mode (rain scan only) during non-critical operations. 


NPOL data in general are considered high-quality. Any issues with data quality will be 
described in the next section. 

3.2.2 Quality Control and Archival of NPOL data 

A tremendous effort has gone (and is ongoing) into the quality control (QC) of the 
NPOL data. Several issues were identified with NPOL data during and after the field 
campaign. Data quality issues include: ground clutter, backlobe and sidelobe echo, wet 
antenna signal loss, and reflectivity calibration. Strong clutter returns occurred during 
strong refractive gradients periods which normally coincided with a minimum in 
convective activity. Strongly contaminated scans were identified by hand and flagged. 
The remaining dataset is being QC’ed for clutter using the standard QC algorithms 
developed by the NASA TRMM office. 

Because of the design of the NPOL antenna, there are occurrences of false echo 
from sidelobe and backlobe returns. A sidelobe signature is seen as a high reflectivity 
circular ring center at the radar location. These scans are easily identified and were 
removed from the QC’ed dataset. Backlobe echo is the result of a small amount of signal 
being transmitted behind the reflector. A return is seen if the echo is 26 dB or higher 
above the minimum detectable signal. The signature of backlobe echo is a mirrored 
image of the strong echo (but 26 dB lower) that is located at the same range but 180° in 
azimuth from the original target. Backlobe echo has been difficult to remove because it 
‘looks’ similar to weak, stratiform precipitation. Currently, we are deploying a 
combination of algorithms that examine the vertical structure and standard deviation of 
the specific differential phase to remove the backlobes. 

Periods when NPOL has a wet antenna cause the signal to be reduced 10-12 dB 
from when the antenna is dry. We have not been able to identify a good QC correction 
algorithm for wet antenna. Therefore, we have removed all scans during wet antenna 
periods (identified by coincident rain gauge data). This removal of data is unfortunate, 
but results in only about 5% loss of the data. 

The remaining data have been analyzed and compared with other datasets to 
determine the relative calibration. A comparison with the surrounding WSR-88D radars, 
the TRMM PR, and internal checks using the polarimetric fields indicated that NPOL had 
a reflectivity bias of -6.5 dB on average (Kucera and Frank 2003; Kucera et al. 2004). 
This reflectivity offset has been applied to the QC’ed NPOL dataset. 

A version 1 QC’ed dataset is available from the UND ftp site: 
lftp.aero.und.edu/public/atmos/crystal face/npol/) . This dataset only has the average 
reflectivity bias correction and a simple clutter rejection algorithm applied. Wet antenna 
and strong sidelobe scans have been removed from the QC dataset. A version 2 dataset is 
currently being processed that will account for the backlobe signatures, a daily 
reflectivity correction, ZDR bias correction, and more robust clutter algorithm correction. 
This dataset (considered final) will be delivered to the C-F archive by the end of 
December 2004. 


Quicklook reflectivity products were created for every scan recorded during the 
C-F campaign. These scans have been archived at UND: 
ftp.aero.und.edu/public/atmos/crvstal face/images/ and at the C-F website in the 
following locations: http://espoarchive.nasa.gov/archive/arcs/crvstalf/images/npol 1/. 
http://espoarchive.nasa.gov/archive/arcs/crystalf/images/npol 2 A 
http://espoarchive.nasa.gov/archive/arcs/crvstalf/images/npol 3£, 
http://espoarchive.nasa.gov/archive/arcs/crvstalf/images/npol 4/. 
http://espoarchive.nasa.gov/archive/arcs/crvstalf/images/npol 5/ . 

3.2.3 Development of Polarimetrically-Tuned Rainfall Products 

In close coordination with co-PI’s Cifelli, Tom Rickenbach (NASA/GSFC), and 
Larry Carey (Texas A&M), we have created polarimetrically-tuned rainfall maps from 
NPOL data. Larry Carey tuned his polarimetric rainfall algorithm for the priority day of 
July 23. During this process, he also created pol-tuned reflectivity-rainfall (Z-R) 
relationships for the four surrounding WSR-88D radars (Key West, Miami, Melbourne, 
and Tampa Bay). Rainfall estimate were created for the five radars. Inclose 
coordination with Tom Rickenbach, rainfall maps from each radar were merged to create 
a composite rainfall map for South Florida. The rainfall products have a spatial 
resolution of 2 km and are mapped to a height of 2 km above ground. Rainfall maps 
were generated every 10 min. Rainfall maps were generated from each radar and merged 
for the remainder of the C-F campaign (July 3-29). These version 1 rainfall products are 
available from the UND ftp site: ftp.aero.und.edu/public/atmos/crvstal face/rainmaps/ . 
Evaluation of the merged rainfall maps with rain gauge data for periods outside of July 
23 shows a rainfall bias between radars. Most of the bias can be associated with the 
NPOL polarimetric rainfall algorithm. The amount of data used for the tuning was 
relatively small, which resulted in a large uncertainty in the rainfall algorithm, especially 
at the higher rainfall rates. Also, because of the remaining QC issues mentioned above, 
there is a larger the normal uncertainty in the rainfall maps. These issues are being 
addressed in the remaining QC work. It is anticipated that a version 2 (final) QC’ed 
rainfall product will be delivered to the C-F archive at the same time or shortly after the 
delivery of the QC’ed NPOL data in December 2004. 

3.2.4 Hydrometeor Signature Studies 

This task is the most complete of the proposed tasks. The main reason is the 
unanticipated extra effort involved in QC’ing NPOL data and the creation of the merged 
C-F radar rainfall products. However, microphysical comparisons have started with the 
initial radar/Citation comparison for the July 16 priority case, which was presented at the 
International Conference on Clouds and Precipitation in Bologna, Italy (Poellot et al. 
2004). This work is currently being expanded to compare bulk microphysical properties 
to the in situ measurements from the Citation for July 16 and other priority cases on July 
1,9, 19, 21, and 23). This work will continue into Fall coincident with the final NPOL 
QC activities. 

3.2.5 Collaborative Studies 

There have been several ongoing studies outside the main tasks proposed in this 
work. There currently is a mesoscale modeling study underway in collaboration with Sue 




van den Heever (CSU) that was initiated from detailed study of observed gust fronts 
during C-F (Frank and Kucera 2003). In an effort to study the evolution of convection 
and associated cirrus anvil, a study is underway to track cloud elements using both NPOL 
and available satellite data using a correlation tracking technique. A preliminary study 
for the July 16 priority day was presented at the Fall 2003 AGU meeting (Kucera et al. 

2003) . Another study was conducted to examine nitrogen generation from lightning 
discharge in thunderstorms observed by NPOL and EDOP during C-F (Ridley et al. 

2004) . Finally, there are ongoing research activities examining the characteristics of 
convection during C-F with NPOL and other independent datasets (Heymsfield and 
Belcher (NASA/GSFC), Tao and Lang (NASA/GSFC), and Zipser and Li (U of Utah)). 
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Table 1. Citation Instrumentation Suite 


Instrument 

CFDC 

"cvT 


ON 

CPI 

HVPS 

FSSP 

PMS 2D-C 


CNC 


TDL H20 


King, Rosemount 


State Parameters 
Winds 


Aircraft Parameters 


Parameter 


Ice nuclei concentrations, filter 
samples 


Ice mass, residual particles 


Mass extinction 


Cloud particle images 


Precipitation particle images 


Cloud droplet sizes 


Cloud particle images 


Condensation nuclei 
concentrations 


Water vapor concentration 


Supercooled liquid water 


Temperature, Pressure, Dew/Frost 
Point 


u, v, w, turbulence 


Position, Orientation, Motion 


Investigator 


DeMott (Colorado State) 
Rogers (NCAR) 


Twohy (Oregon State) 


Gerber (Gerber Scientific) 


Heymsfield (NCAR) 


Heymsfield (NCAR) 


Poellot 


Poellot 


Poellot 


Poellot 


Poellot 


Poellot 


Poellot 


Poellot 


Table 2. Citation Mission Summar 



Hours Comments 


July 3 


July 7 


July 9 


July 11 


July 1 1 


July 19 


July 21 


July 23 


July 25 


171444 


175842 


175517 


142548 


180651 


183744 


165043 


183648 


183500 


185225 


150449 


202909 


193833 


171101 


211650 


214245 


214935 


165155 


221418 


225329 


203043 


224621 


224852 


232229 


193147 


235643 

220733 

212418 



Steps and spirals 

Steps and spirals; dense cirrus 


Steps over east and west sites. Mixed 
phase over west site 


Altostratus, Terra overpass 


Anvil steps 


Spirals, steps in anvils 


Ice nucleation, anvil spiral 


Anvil steps, spiral 


Anvil steps, spiral 


Anvil steps, spirals 


Steps and spirals; samples over sites 


Melting layer spiral descents 


Growing tower samples; dust on descent 


Anvil steps, spiral; dust on descent 













































































Appendix A 
Mission Summaries 



Project Crystal-Face 
July 3, 2002 - July 29, 2002 
Citation Flight Summaries 


Notes 


Segment: A section of the flight where relevant samples where taken. 

Three different types of segments: 

1) Sample Line 

2) Spiral (up or down) 

3) Step Climbs (ascending or descending) 

Leg: The level sections of a step climb where relevant samples where 

taken. 

Cloud Region: The cloud or area of clouds where the sample was taken 

from. 

Cloud regions are numbered and segments with 
same cloud region numbers indicate that those segments 
are in the same cloud region. If NC is shown in the 
cloud region column it indicates that there were no 
clouds present. 

Parameters: Six different parameters included in the summaries: 

1) Temperature given in degrees Celsius 

2) Frost Point given in degrees Celsius 

3) Wind Direction given in degrees 

4) Wind Speed given in meters per second 

5) 2-DC Concentration given in number per liter 

6) Maximum 2-DC Concentration given in number per 
liter 

7) 2-DC Median Volume Diameter given in microns 

Values for each parameter are averaged over the segment 
or leg if the segment type is a step climb. 

The temperature and frost point values are given as a 
range of the lowest to highest temperature for spiral 
types of segments. 

The frost point parameter is measured by the laser 
hygrometer, unless other wise noted. 

During a spiral type segment there is no valid wind 


speed or direction. 


Additional notes pertaining to specific flights are given on the spreadsheets. 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/11/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/16/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/18/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/19/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/21/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/23/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averged over segment time.) 

Date: 7/25/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/26/02 
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Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/28/02 
































































Flight Summary 

(Flights broken up into segments, step climbs broken down further into individual legs) 
(All parameters are averaged over segment time.) 

Date: 7/29/02 
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